Adiabaticity Conditions for Volatility Smile in Black-Scholes Pricing Model 
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Our derivation of the distribution function for future returns is based on the risk neutral approach 
which gives a functional dependence for the European call (put) option price, C(K), given the strike 
price, K, and the distribution function of the returns. We derive this distribution function using for 
C(K) a Black-Scholes (BS) expression with volatility, a, in the form of a volatility smile. We show 
that this approach based on a volatility smile leads to relative minima for the distribution function 
("bad" probabilities) never observed in real data and, in the worst cases, negative probabilities. We 
show that these undesirable effects can be eliminated by requiring "adiabatic" conditions on the 
volatility smile. 
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I. INTRODUCTION 

One of the simplest "products" on the derivative fi- 
nancial market is the European call (put) option [HE]- 
Considering the risk neutral approach, the price of the 
European call option, C = C(St, K, T, r), is defined by 
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-tT 
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(1) 



where St is the stock price at time t = T, K is the strike 
price of the option, T is the expiration time (time to ma- 
turity) of the option, r is the interest rate and P(St) > 
is the distribution function of the stock prices in a "risk- 
neutral world" (J °° P(S T )dS T = 1). 

Eq. ([IJ is too general since it does not place any re- 
strictions on the underlying stock price distribution func- 
tion, P(St)- To calculate explicitly the option price, C, 
using Eq. (fTl), one must know the distribution function, 
P(St)- Consequently, one must make some assumptions 
about the stock prices. An important achievement in the 
theory of option pricing is the Black-Scholes (BS) theory 
which gives analytic solutions for the European call and 
put options [3]. 

In particular, for the European call option, a solution 
of the BS equation is given by Eq. ([IJ , if one assumes for 
the distribution function, P(St), a log-normal distribu- 
tion, 
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where x = In (St / S (t)) — r(T —t) is the logarithmic return 
deprived of the risk-free component, S(t) is the stock 
price at time t and a is the stock price volatility. For 
seek of simplicity, in the following we consider t = and 
we define S = S(t = 0). Substituting Eq. (JsJ) in 
an explicit expression for the price of the European call 



option which satisfies the BS equation [3J is obtained, 

C BS = SoNfa) - Ke- rT N(d 2 ), (3) 

where 
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The distribution function ^ follows from a stochastic 
model for stock prices, 



dS = rSdt + aSdz, 



(5) 



where dz is a Wiener increment [4]. It can be shown it 
is never optimal to exercise an American call option on a 
non-dividend-paying stock early pQ , [5] ; therefore Eq. ^ 
can also be used to estimate the fair value for this kind 
of options. 

There are some problems with the expressions for C 
given by Eqs. Q-([3|. Indeed, one can derive any option 
price from Eq. (IT]) , using different assumptions about the 
distribution function, P(St)- To derive from Eq. (JlJ a 
result for C which will even approximately coincide with 
the real market price, C M , one must specify a distri- 
bution function for future stock prices, P(St)- On the 
other hand, the expression given by ([3| is (a) too specific, 
and (b) derived using rather strong restrictions. Namely, 
the stochastic process Eq. ^ does not account for cor- 
relations of returns, x and, moreover, the volatility, a, 
and the interest rate, r, are not well-defined parame- 
ters (given the actual data). As a result, the expression, 
C BS , often does not coincide (even approximately) with 
the corresponding market option price, C M . Useful ap- 
proaches have been developed which partially solve the 
problems mentioned above. 
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We shall mention here one analysis which is related to 
that presented below. This analysis deals with building 
"implied trees" [BJ. There are many variations of this 
approach, but the main idea is based on the solution of 
the inverse problem: a search for a stock price model that 
corresponds to the real market prices of options, C M . A 
more restricted problem is to search for a stock price 
model that effectively deals with the volatility smile. 

In this case, one starts with the BS formula ([3|, (even 
for American options) but instead of choosing a fixed 
volatility, a = constant, one uses the dependence, a = 
o~(K) (volatility smile). To some extent, this dependence 
corresponds to the "real behavior" of the volatility, a, if 
one wants to use Eq. (j3| as the "zeroth approximation" 
for option pricing. Details for building trees (including 
implied trees) for stock prices can be found in [7HTU]. 

There are still some problems with these trees. For 
example, the corresponding "implied" stock prices, St, 
can have "bad" (negative) probabilities which must be 
eliminated. A solution to this kind of problem, despite 
the simplicity of the calibration procedure, was proposed 
in Refs HU[T2]. 

In this paper, we discuss the inverse problem for Eq. ([!]) 
using the following approach. First, using Eq. ([!]), we 
build the distribution function for future stocks prices 
and returns from the empirical data for the market op- 
tion prices, C(K). Second, we build the returns distri- 
bution functions using for C(K) a BS expression, C BS , 
with a volatility, a — <j{K), in the form of the volatility 
smile. In particular, we show that the condition of the 
absence for relative minima in the probability distribu- 
tion function (PDF) of returns, (or elimination of "bad" 
probabilities) leads to the condition of "adiabaticity" for 
the volatility smile. This condition can be introduced 
in the fitting procedure of the volatility smile to get a 
probability returns distribution more similar to the ac- 
tual one. In this way one can avoid to generate arbitrage 
opportunities (negative probabilities) in the option pric- 
ing methodology (exotic derivatives) and can get a more 
reliable estimation of the implied volatility. The latter 
has a key role in the scenario generation and in value 
at risk (VaR) estimation and has application in the risk 
management activities. 



II. THE INVERSE PROBLEM FOR THE 
STOCK PRICE DISTRIBUTION FUNCTION 



In this Section, we derive an explicit expression for the 
distribution function for the future stock prices, P{St) 
and for returns P(x). In Eq. ([I]) the distribution func- 
tion, P(St), can be rather arbitrary but it is natural to 
assume that P(St) does not depend on the strike price, 
K. According to Eq. ([I]), the option price, C, is expressed 
explicitly through the strike price, K. Differentiating C 
in Eq. twice with respect to K, we have [13) . 
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In Eq. d6j), we indicate only the dependence C(K) in 
the option prices. In particular, applying Eq. ^ to 
C BS given in Eq. ^ we derive a distribution function, 
P B (St), which we present in the form, 
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Analogously, the distribution of returns for the Black- 
Scholes model is Gaussian, as expected: 
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where x Q — a 2 /2T. 

We can try to consider the inverse problem substituting 
the dependence, a = o(K), in Eq. ^ and evaluating the 
distribution of future stocks price and returns, applying 
After differentiation we get: 
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Obviously, we can get the expression for the distribution 
of returns by a simple change of variable 
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(11) 



K=Si 



so that, 
P*(x) 



1 



V2ira 2 T 



exp 



(x + x ) 2 
2a 2 T 



F(x;T,a), (12) 
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FIG. 1: a) Volatility smile given by Eq. (171, with the fol- 
lowing parameters: g = 0.1, T = 0.5, n — 0.04, \ — 2.7. b) 
Log-distribution given by Eq. (|12[) for the volatility smile a). 



where, in a similar way, we have defined: 
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From ( 13 ) it is clear that if a is constant ^ and (12]) are 
the distributions for the standard Black-Scholes model; 
we will call them "zeroth-approximation" distributions. 
If a 7^ const, the term F(x;r,T,a) could "perturb" the 
relative zero approximation (Gaussian) giving rise to dis- 
tributions that cannot fit real data. As we will show here, 
it is possible to get distributions with relative minima 
(not observed in real returns distributions) and, in the 
worst case, negative probability. 



The rest of the paper is organized as follows. In Section 
|III[ we analyze FX market data and we propose a suit- 
able function to fit the volatility smile. We also determine 
the typical range of the parameters we use to perform the 
fit. We will use these ranges to check our adiabatic con- 
dition. In Section IV we show qualitatively the reason 



why we have the presence of the relative minima in the 
returns distribution and why an adiabatic approach can 
be suitable to describe the problem and used to avoid 
these "bad" probabilities. In Section |VJ we show our 
numerical and theoretical results about the adiabaticity 
parameter. We also show that there is a critical value of 
the adiabatic parameter which can be used to determine 
whether the returns PDF will have relative minima. We 
discuss the relation of this critical value to the other pa- 
rameters of the fit. Finally, in Section [Vl| we present our 
conclusions. 



III. THE VOLATILITY SMILE: REAL MARKET 
DATA 

Typically, traders on option markets and practition- 
ers consider the BS model as a zeroth order approxima- 
tion that takes into account the main features of options 
prices. To get a pricing closer to the actual data, they 
consider the volatility as a parameter that can be ad- 
justed considering the inverse problem given by Eq. ^ 
and the real price of call and put options. In this way 
a more reliable value of the volatility (implied volatil- 
ity) can be obtained and it can be used to price more 
complex options for which analytical solutions are not 
available. The value of the implied volatility depends on 
the value of the strike, K, in a well-known characteristic 
curve called the smile volatility (typically for foreign cur- 
rency options) whose shape is approximately parabolic 
and symmetric, or skew volatility (typically for equity 
options) when asymmetric effects dominate [TH - H7] . 

An intuitive explanation of this shape can be found 
if an actual returns distribution is considered. In fact, 
it is well known that the tails of the returns PDF 
are not Gaussian but exhibit a power law decay (fat- 
tails) [T5J[TS]. On the contrary, BS model assumes that 
the PDF of returns is Gaussian thus underestimating the 
actual probability of rare events. To compensate for this 
model deficiency, one has to consider the greater implied 
volatility for strike out of the money then for strike at 
the money. 

In this paper we focus our attention on the volatil- 
ity smile of foreign currency options and we neglect the 
skew effect |20) . To perform our analysis we consider 
the volatility smile as a function of the A of the option 
(defined by Eq. (14)), the time to maturity, T, and the 



currency considered. We consider specific days, for which 
volatility is not affected by the skew effect, and we use 
Bloomberg as data provider. In the BS model, the A of 
a call options is defined as: 
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Inverting this relation is possible to get an expression for 
x: 

x = o- 2 /2T-o-VT$- 1 (A), (15) 
where ^~ 1 (x) is the inverse of the error function: 
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In Fig. [2] we show an example of volatility smile in terms 
of our variables and a suitable fit given by the function: 



<t(x) = g 



1 + (X - 1) 



(x + g 2 T/2) 2 
(x + g 2 T/2) 2 + n 



(17) 



where g, x, n are fitting parameters. In this case, g rep- 
resents the minimum of the volatility smile, \fri is the 
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FIG. 2: Typical volatility smile and the relative fit obtained 
with (171. The parameters of the fit are: g = 0.1758(5), 
X = 1.20(9), n = 0.00030(9). We get the real data using 
Bloomberg provider and they refer to the AUDUSD currency 
with time to maturity T = 1/365 years. 



half width at the half height, while g(x — 1) represents 
the height of the smile. In particular x is the ratio be- 
tween the limiting value of a and g as x approches oo. In 
this way the variation of a is bounded between g and g\- 
In the light of the intuitive explanation of the volatility 
smile proposed above and since from ^ it follows that 
the average value, 



a 2 T 
2 ' 



(18) 



one expects that the minimum of the implied volatility 
occur at x = —g 2 T/2 as required by our fitting function. 

Repeating many times the interpolation procedure 
considering different values for T and currencies, we can 
determine typical parameters that can fit a wide range 
of volatility smiles; in the following we will use this in- 
formation to check our results. 

Let us notice that the following relation between n, g 
and T holds: 



n oc Tg 2 



(19) 



as shown in Fig. [3j This gives a scaling rule that can be 
used to determine the range of n, fixing T and g. 



IV. FIRST APPROXIMATION OF THE 
VOLATILITY SMILE: THE SQUARED WELL 

In this Section we show qualitatively the reason why 
there is a relative minimum in the returns distribution 
and why an adiabatic approach can describe the prob- 
lem of avoiding these "bad" probabilities. To keep the 
problem simple, we consider, as a first step, a volatility 
smile modeled by a squared well defined as follows: 



a{x) 



a i for \x\ < X\ 

a 2 otherwise 



(20) 



FIG. 3: Relation between the parameters n, g, T. We fit 
72 volatility symmetrical smiles (Bloomberg) considering dif- 
ferent currency (EURUSD, AUDUSD, EURCHF, EURGBP, 
EURJPY, GBPJPY, GBPUSD, USDCAD, USDCHF, USD- 
JPY) and time to maturity (1 day, 1-3 weeks, 1, 2, 3, 4, 6, 
9 months, 1, 1.5, 2, 3, 4, 5 years) with the function (17 1. 
We dowload the data on 21/10/2009 and on 01/02/2010. 
We also show the best linear fit ln(<; 2 T) = ln(n) + c, where 
c = -1.95(12). 
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FIG. 4: a,b) Discontinuous squared well as a fit of the volatil- 
ity smile. c,d) The corresponding PDF with (c) or without 
(d) spurious minima. 



where <T2 > oi and x\ are positive constants. See 
Fig. |4]a,b. The distribut ion functions corresponding to 
two values of <7j, i = 1,2 are shown in Fig. |4];,d. Indi- 
cating as ±xf the abscissa of the intersections between 
the two distributions, it is clear that a sufficient condi- 
tion for avoiding spurious minima is xi < x\. A rough 
estimation of x\, ignoring the term, xq, usually small, is 




(21) 



where \ = fa/o'i- Therefore, a sufficient condition to 
avoid minima in the PDF is to use, as a fitting function a 
square well depending on the parameters X\,a\,a2, such 
that X\ < x\. Therefore, a standard fitting procedure of 
the volatility smile with a square well, constrained by the 
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minima, while this does not happen in the region to the 
right of it. It is then clear that for a given set of fixed 
parameters (<?,T), there is a relation between x an d n 
that allow one to obtain a PDF without minima (min- 
ima are not observed in real data). The main goal of this 
paper is to determine a simple relation that determines 
whether or not the parameters of a volatility smile fit are 
consistent with real returns distribution and if they could 
give a reliable option pricing. 
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FIG. 5: a) Critical ratio \ c = 02 /&i as a function of n for 
fixed T = 0.5, and different g as indicated in the legend. The 
arrow indicates the direction of growing g. In the region to 
the right of the lines the PDF do not have minima, while 
in the left hand region it has. b) The same as a) but for 
fixed g = 0.1 and different T values. The arrow indicates the 
direction of growing T. 



V. NUMERICAL AND THEORETICAL 
RESULTS 

In this Section we show our numerical and theoretical 
results about the relation between the set of parameters 
n,g,T and the critical adiabatic parameter Xc- Using a 
numerical simulation, wc kept fixed n,g,T and we con- 
tinuously increased the parameter x until we found a 
zero-derivative point for some x ^ — g 2 T/2. In this way 
we could determine numerically the critical Xc- We re- 
peated this approach for a wide range of the parameters 
values, as shown in Table [TJ where we used the parameter 
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condition x\ < x\ solves the problem of avoiding spurious 
minima, even if is very rough. 

If we consider a volatility smile with a continuous vari- 
ation from &\ to (72 1 we can get, instead of a discontinu- 
ity, the onset of a relative minimum. The latter can be 
avoided if the variation between a\ and a 2 is slow enough 
so that the connection between the two PDF's takes place 
keeping constant the sign of the first derivative of the dis- 
tribution during the whole transition. So there will be a 
critical "speed" of the transition that will generate zero 
derivative points which will not correspond to the maxi- 
mum of the distribution. In this case the variable related 
to the time is Xx, while x — 02/01 can be identified as a 
"distance" . To be more precise, one should consider that 
for x < x\, P ai {x) > P a2 (x), so that the effective "time" 
should be: x\ — x\. 

We will define an adiabatic transition as one whose pas- 
sage from o\ to <72 is sufficiently slow, so that relative 
minima in the distribution are not generated. In this 
same spirit, we define the critical adiabatic parameter 
% c (n, </,T) as the minimum value of x that generates a 
minimum in the distribution, Eq. (12). This effect has 
been shown in Fig. ([5]), modeling the volatility smile us- 
ing (17), so that (Ti — g, 02 = gx- We fix the parameters 



g, T and we vary n and x, seeking relative minima in the 
PDF. The lines in Fig. [5] obtained respectively at fixed 
T (a) and fixed g (b), divide the plane of parameters into 
two regions: to the left of the lines the PDF has spurious 



TABLE I: Range of the parameters of the numerical simula- 
tions 



p = n/g 2 T instead of n, due to the scaling relation ( 19 ). 
In order to obtain the relation Xc — /t(^,5), we use the 
following fit function: 



fr(n,g) = a 



9 2 T 



(22) 



This has been obtained assuming that the value of the 
critical parameter Xc depends on the rescaled "time" of 
the transition (in our model given by p). We also consider 
a further term j^/Tgp^ to take into account the time 
correction, x c , as explained in Sec. (IV). In this case we 



make explicit the dependence of the time correction on 
T and g as suggested by a 2 in Eq. pTJ. In Eq. p2l, 



a, /?, 7, 5 are the fitting parameters whose values are given 
in Table HU 

In Fig. [6] we show the result of our fit for a few selected 
values of g and T. 

The whole procedure can thus be summarized as fol- 
lows : 

• The real volatility smile, usually given for a fixed 
T can be fit by a function dependent on three pa- 
rameters g, n, X-, as indicated in Eq. (17) and the 
optimal values g opt ,n opt , 



X° pt are returned 
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1.4373 ± 0.0002 


13 


0.2787 ± 0.0006 
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-0.1738 ±0.0002 
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0.4683 ± 0.0006 


mean squared errors 


1 x 10" 5 



TABLE II: Fitting parameters and relative errors. 
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a suitable generalization to include the volatility smile 
effect. Then we considered the inverse problem and we 
analyzed the relative returns distribution, Eq. ( |i"2| , vary- 
ing the typical parameters of the volatility smile. We 
showed that, for some values of the parameters, it is pos- 
sible to get relative minima in the returns distribution 
(bad distribution) that are never observed in real distri- 
butions. We demonstrated that bad distributions can be 
eliminated by requiring adiabatic constraints (intuitively 
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FIG. 6: Critical adiabatic parameter as a function of n for few 
selected pairs of values of (g, T) as indicated in the legend. 
The points are numerical data, straight lines are the result of 
fitting procedure. 



FIG. 7: a) Volatility smile as a function the returns x. Dots 
indicate real data, red curve is the non adiabatic fit, while the 
blu one represents the adiabatic (constrained) fit. b) PDF of 
returns for the two curves indicated in a). 



The optimal values 



ical Xc 



gOpt^Opt 



are inserted in 



Eq. (22 1, with a, j3, 7, S given in Table II and a crit- 



obtained. 



If x° pt < Xc then we know that relative minima 
in the PDF do not exist. Otherwise we should per- 
form a fitting procedure to the volatility smile using 
Eq. (17), constrained by x < Xc- 



An example of the previous procedure has been shown 
in Fig. [7] where the PDF with unwanted minima and the 
"corrected " one is shown together with the correspond- 
ing fitting curve to the volatility smile. As one can see 
the price to pay in order to get a smooth PDF is very 
small: the two fitting curves for the real volatility smile 
are similar, but the PDF has, in the latter case a more 
realistic behavior. 



VI. CONCLUSIONS 

We started from the pricing equation of the Black- 
Scholes model for an European call and we considered 



justified with the example of the squared well) on the 
volatility smile and we gave a numerical formula to de- 
termine the value of the adiabatic critical parameter, Xc- 
In this way we provide an easy-to-use tool to determine if 
a volatility smile fit is consistent with the model hypoth- 
esis (P(x) > 0) and if it can generate a suitable returns 
distribution. A reliable estimate of the implied volatility 
has application in the risk management activities and in 
the pricing of exotic derivatives, where, in general, the 
implied volatility is an input of more complex models. 
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